ORNL-CTP-9805 and 



[iep-ph/9809497 



Drell-Yan Cross Section and J/ij] Production in High-Energy 
Nucleus-Nucleus Collisions 



, Cheuk-Yin Wong 

^\ • Phystcs Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA 

0^ ' E-mail: wongc@ornl.gov 



We use the Drell-Yan differential distribution dN^y /dEx in high-energy nucleus- 
nucleus collisions to obtain a relation between the transverse energy Et and the 
impact parameter b. Such a relation is then utilized to study the transverse-energy 
dependence of J/ip production in Pb-Pb collisions, using the J/ip absorption model 
presented previously. The anomalous Pb-Pb suppression data at 158A GeV can be 
explained if one assumes the occurrence of a new phase of strong J/ip absorption 
' when the energy density exceeds 4.2 GeV/fm^. The results are extended to make 

^ ' predictions for J/^p production at higher collision energies. It is found that J/i/> 

, survival probabilities reach the lowest survival limit when the nucleon-nucleon 

■ center-of-mass energies is greater than about 35 GeV. 

i 1 Introduction 

OO ■ n 

0^ ' Following the initial suggestion by Matsui and Satz 111 that J/^p production, i, 



X 



11 



suppressed in a quark-gluon plasma, the NA50 experimental observation!: 
pi i' of an anomalous J/j!/i— absorption in Pb-Pb collisions has led to a flurry of 

' . theoretical activity eHIj. The absorption is anomalous in the systematics of 

^j' the total normalized yields as a function of the sum of the radii of the colliding 

(~| , nuclei. It is also anomalous in the differential J/ip survival probability, which 

shows a discontinuity as a function of the transverse energy. A central question 
is whether these discontinuities indicate the occurrence of a new phase of strong 
J/ijj absorption, as expected when a quark gluon plasma is produced. 

Theoretical studies on J/^ suppression require the relation between the 
transverse energy Et and the impact parameter 5 of a collision. We shall first 
show how we can obtain such a relationship from the Drell-Yap distribution. 
We shall use the J/V' absorption model presented previouslylla to study J/ip 
absorption in Pb-Pb collisions at piab = 158A GeV (^s = 17.3 GeV). The 
results are then utilized to predict J/ ip production at higher energies where we 
find interesting trends as the collision energy increases. 

2 Drell-Yan Cross Section 

In NA50 experiments, the transverse energy is measured in coincidence with 
the detection of Drell-Yan and J/tp dimuon pairs. Based on the assumption of 
a monotonic dependence, the relation between Et and b can be easily obtained 
from the Drell-Yan differential cross section da^y /dET or the Drell-Yan dis- 
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tribution dN^^/dEr which is equal to {da^^ / dET)[N^§-{totoJlln^^ {total)]. 
One relies here on the experimental nuclear mass dependencesE3i3 which indi- 
cate that initial- and final-state interactions have only very small efi^ects on the 
Drell-Yan cross section in nuclear collisions. Neglecting initial- and final-state 
interactions, the Drell-Yan cross section element in the collision of nucleus A 
with nucleus B is given by 

da^^^ABTAB{b)dba^jPy, (1) 

where TAsib) is the thickness function, 

TAB{b)^ J dbATA{bA)TBib^bA), (2) 

and TAibiX and Tsibs) are the thickness functions for nucleus A and B re- 
spectivelyim. Therefore, we have 

if = AB^, T^sm ^. (3) 

From these relations, we can get a relation between Et and b given by 



(4) 



where Y.{Et ^ oo) = I](6 = 0) = 1 and 



/•oo J AB 

o^^itotal) = ^dE'^. (5) 

If the Drell-Yan distribution dN^y/dEx is measured instead, the function 
T,{Et) can be expressed in terms of dN^y / dEr by 



^(Et) ^ -ABk^ I -T^dE'r (6) 



rEr ^^A B 

N^^ (total) Jo ~dE'j. 

where j-ca jtsjAB 

N^^itotal) = ^dE'^- (7) 

To obtain Et as a function of fo, one can use a graphical method where one 
obtains 'S{Et) from experimental data and from TabI^) which can be 
calculated from the known geometry of the colliding nuclei. One then uses h 
and Et as the abscissa, and plots S as the ordinate, resulting in the curves 



2 



of 'S{Et) and From Eq. the horizontal line of constant S then 

intercepts the two curves at b and its corresponding Exib). Alternatively, one 
can parameterize Exib) as 

Erib) = — 1 , (8) 



1 + exp[(6 - bo)/a] ' 

which gives db/dEx. One then uses Eq. (^) to determine da^y/dET- The 



Drell-Yan distribution for Pb-Pb collisions calculated with Etq 
bo = 6.7 fm, and a = 3 fm is shown as the dashed curve in Fig. 1. 

160.0 



135 GeV, 



I 



140.0 
120.0 
100.0 
80.0 
60.0 
40.0 
20.0 



O Pb+Pb (NA50 1996) 

- - No Ej spreading 

— 'With Ej spreading 




150 



(GeV) 



Fig. 1. The Drell-Yan distribution as a function of the trans- 
verse energy in Pb-Pb collisionSpaL.158A GeV. The data points 
are preliminary data from NASOcrEa. 

One notes that for a given value of 6, the Et value cannot be a sharp 
distribution. One can consider the righthand side of Eq. (H) as defining the 
mean value £'t(6), and include the spreading of Et about this mean value by 
introducing a distribution function D(b,ET), 

D{b, Et) = ^-[Er~EHb)?m-m\ (9) 

\/27rcr(6) 

Then, with the inclusion of the spreading of Et, da^y/dET is given by 

^ = ABaPPy J db TABib) Dib, Et). (10) 
From the results of Ramello et alB'El, a{b) in Eq. (^ can be parameterized as 



a{b) 



1 + exp[(fo - b„)/aa] 



(11) 
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where ctq = H GeV, 5^ = 10 fm, and a^r = 2 fm. The Droll- Yan distribution 
with the inclusion of Et spreading is shown E^^riiic solid curve in Fig. 1 , which 
agrees well with the preliminary NA50 dataElEa. The Drell-Yan distribution 
provides a relation between Et and b which is needed to analyze the transverse- 
energy dependence of J/ip and ijj' production. 

3 Model of J/ip and tp' Absorption 

We can use the J/ip absorption model presented previously in [111] to study 
J/ip and ip' absorption. In this model, a produced J/ip (or its precursor) meets 
the projectile and target nucleon at high energies and is absorbed with a cross 
section (Jahs[J li^—N). It also collides with produced hadrons (comovers) at low 
relative energies and is absorbed with a cross section aahs{J/ip—h). In a nucleus- 
nucleus collision, each nucleon-nucleon collision is a possible source of J/ip and 
■0' precursors. It is also the source of a fireball of produced hadrons which can 
absorb J/ip and ip' precursors produced by other nucleon-nucleon collisions. 
One follows the space-time trajectories of precursors, baryons, and the centers 
of the fireballs of produced hadrons. Absorption occurs when the space-time 
trajectories of the precursors cross those of other particles. Using a row-on- 
row picture in the center-of-mass system and assuming straight-line space-time 
trajectories, we^obtain the differential cross section for J/iJj production in an 
AB colhsion asEj 



rf%^(b) f dh. 



1 - 



i-rjb-bjaabs(J/^-iv) 



F(b^,b), (12) 



where F{hA,h) is the survival probability due to soft particle collisions. To 
calculate F{hA, b), we sample the target tranverse coordinate hA for a fixed im- 
pact parameter b in a row with the nucleon-nucleon inelastic cross section cri„ . 
In this row, BTsih — byi)cri„ projectile nucleons will collide with ATA{hA)(Jin 
target nucleons. We construct the space-time trajectories of these nucleons to 
locate the position of their nucleon-nucleon collisions. These collisions are the 
sources of J/tp and ip' precursors and the origins of the fireballs of produced 
particles. For each precursor source from the collision j and each absorbing 
fireball from the collision i at the same spatial location, we determine the time 
i'^- when the precursor source coexists with the absorbers in the state of pro- 
duced hadrons. The survival probability due to this combination of precursor 
source and absorber is then exp{— fc^/ji^}, where the rate constant k^h is 

ki^h = Ph^VhCTiihsiJ/^p-h), (13) 
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Vfi is the average {J/ij^yh relative velocity, p^'^ is the average produced hadron 
number density per NN collision, 

NN _ "^^/i ^ /I.N 

dy aUd/lY ^ ' 

dN^^ / dy is the particle multiplicity per unit of rapidity at j/^j^j = for an 
NN collision, d ^ 2.46 fm is the internucleon spacing, and 7 = •\/s/2mnucicon 
is the LorentZpContraction factor. We can use the parametrization of Eq. (2.4) 
in Reference [113] to obtain dN^^ /dy as a function of the collision energy. 

When we include all possible precursor sources and absorbers, F(b^,b) 
becomes w< , , n n 

n=l ^ ^ j=l i=\^i=^j 

where A^>(b^) and iV<(byi) are the greater and the lesser of the (rounded- 
off) nucleon numbers ATA(byi)cri„ and BTsih — b^)cri„, a{n) = 2 for n = 
1, 2, 7V< - 1, and a(A^<) = A^> - iV< + 1. The function 6 is zero if A = 1 or 
B ^ 1 and is 1 otherwise. 

We assume that the Bjorken-type longitudinal expansion occurs after the 
last nucleon- nucleon collision takes place at that spatial point. The survival 
probability F can be determined from k^h and 

^tr,+tf- MaxiU + th, tj + t,g). (16) 

Here we shall set the cc production time tec equal to 0.06 fm/c, the hadroniza- 
tion time equal to 1.2 fm/c, and tf is related to the freezeout time tfrcezcout 

by 

tf = t,, + tn ln(ifreezcout/i/.) - tn + t^ Hp^,"" / p^Z..out)- (17) 

We use a freezeout density of Pfroracout = 0-^ hadrons/fm^, corresponding to 
a produced hadron freezeout separation of 1.3 fm. The if)' production cross 
section can be obtained from the above equations by changing J /ip into V''- 

We have used this absorption model tn stnHy the experimental J /ip and ijj' 
data in p-A and nucleus-nucleus collisions QEj. We found that J/ip absorption 
by produced hadrons, as revealed by 0-Cu, 0-U, and S-U collisions, is small, 
and the experimental Pb-Pb yield is much smaller than the extrapolated results 
if one assumes J absorption by nucleons and produced hadrons only. 

The deviation of the J/ip data in Pb-Pb collisions from the conventional 
theoretical extrapolations of p-A, 0-A, and S-U collisions suggests that there 
is a transition to a new phase of strong absorption which sets in when the 
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local energy density exceeds a certain threshold. We have extended the ab- 
sorption model to describe this traDsition in terms of the critical number of 
NN collisions at a spatial point EnlJ. We can reformulate the model here in 
terms of the critical energy density. Evaluated in the nucleon-nucleon center- 
of-mass system, the energy density at the spatial point, which has had v prior 
nucleon-nucleon collisions, is approximately 



NN 

e = vp,^ mt 

where mt is the average transverse mass of a produced hadron. 
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Fig. 2. The J/ip distribution as a function of the transverse 
energy in Pb-Pb collisions a±_158A GeV. The data points are 
preliminary data from NASOQ'EJ. 



We postulate that soft particles make a transition to a new phase with 
stronger J /ip absorption at a spatial point if the energy density at that point 
exceeds ecrit- When this happens, the matter in the new phase absorbs J /ip 
with a greater strength. To keep track of the same matter density, we still 
use the same number of particles in the new phase as in the hadron matter, 
but the particles in the new phase absorb J /ip with an effective cross section 



CahsiJ/ip-x)- The quantity k^htfj in Eq. ( [T^ ) becomes k 
the new rate constant k^^x is 



Ph^vx<JiihsiJ/i>- 



k^xtfj, where 



(19) 



and the quantity tfj is the time for a J/ip produced in collision j to coexist at 
the same spatial location with the absorbing soft particles produced in collision 
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i in the form of the new phase. The duration when the matter is in the new 
phase can be inferred by assuming that the matter undergoes Bjorken-type 
longitudinal expansion after all the nucleon-nucleon collisions have occurred 
at that point. Because the soft particles are relativistic, we take all relevant 
velocities v in Eqs. ( p^ and (|l^) to be 1. 

We obtain good agreement with Pb-Pb data using the following parame- 
ters: aahsi'fp-N) = aabs{ip'-N) = 6.36 mb, e^it = 4.2 GeV/fm^, crabs(J/V'- 
h) = 0.15 mb, aa.hs{J/'^~x) ~ 3.7 mb, and aahs{ip' ~h) = 11 mb. Note that 
(Jahs{J/i^—x) » anhsiJ/ip—h). We show in Fig. 2 the J/V' distribution in 
Pb-Pb collisions at 158 A GeV. The solid curve is the result when one assumes 
a new phase and Et spreading, the dashed curve is the result with the new 
phase but without Et spreading. The dashed-dot curve is the result without 
the new phase, but the Et spreading is included. Ther£ is good agreement of 
the theoretical curve with the preliminary NA50 datau when one assumes a 
new phase of strong J/ip absorption and takes into account the Et spreading. 




Ej. (GeV) 

Fig. 3. The ratio of Ba{J/'ip)/a{DY) for Pb-Pb collisions at 
158 A GeV as a function of Et- The data points are preliminary 
data from NA5ol. 

In Fig. 3 we show Ba^^{J/^) ja^^ (DY) for Pb-Pb collisions at 158A GeV. 
The short-dashed curve gives the results when there is absorption by nucleons 
only, while the long-dashed curve gives the results when one assumes absorption 
by nucleons and produced hadrons. Depending on the critical energy density, 
the ratio Ba'^^{J/ip)/cr^^{DY) deviates from the long-dashed curve (for no 
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new phase) at different transverse energies. For tcrit = 4.2 GcV/fm"^, the break 
occurs at Et ^ 40 GcV which is close to the break observed experimentalljd. 
For Ecrit — 3.0 GeV/fm^, it occurs at Et 20 GeV, and for edit = 5.5 
GeV/fm'^, it occurs at Et ~ 90 GeV. However, the discontinuity in the the- 
oretical resiilts is much smoother than the discontinuity in the preliminary 
NA50 datala shown in Fig. 3. 

4 Dependence of J/tp Survival Probability on Collision Energy 

The above formulation allows one to calculate the J/ip survival probabilities, 
represented by S{b) = [Ba^^{J/^)/a^'^{DY)]/[BaPP{J/ip)/aPP{DY)], as a 
function of the impact parameter for different collision energies. The results 
are shown in Fig. 4. The short-dashed curve is the survival probability when 
J/ip is absorbed by nucleons only. The dash-dot and the solid curves are 
the results when there is a new phase. The dashed curves are the survival 
probabilities when there is no new phase. The numbers label the collision 
energy, ^s. 
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Fig. 4. Theia.tio [Ba'^^{J/^^J)/a'^^{DY)]/[BaPP{J/^^:)/aPP{DY)] 
for Pb-Pb collisions at different collision energies as a function of the 
impact parameter. 
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Figure 4 shows that if a new phase of strong J /ip absorption is assumed, 
the J /ip survival probabihty reaches the "lowest survival" (LS) limit (the solid 
curve) when the nucleon-nucleon center-of-mass energy ^/s reaches about 35 
GeV. This is the limit at which the energy density of the produced matter is 
so large that no J particles, except those produced at the peripheral region 
of the collision, can survive the nucleus-nucleus collision. If there is no phase 
transition, the lowest survival limit is reached at about ^/s ^ 70 GeV. 

The existence of this lowest J /ip survival limit shows up as a universal 
curve when one plots the J / ip survival probability as a function of b for different 
energies. This can be obtained by transforming S{Et) to S{b), utilizing the 
method of Section 2 to make a correspondence between Et and b by means 
of the Drell-Yan distribution. Upon verification of the existence of the lowest 
survival limit, the energy at which this LS limit is reached then gives a good 
indication as to whether the new phase of strong absorption occurs or not. 

The above results have an important implication for the RHIC collider. 
At the energy y/s = 200 GeV, one expects that the J/ip survival probability, 
as measured by S{b) = [B(T^^{J/ip)/a^^{DY)]/[BaPP{J/ip)/aPP{DY% will 
be in the lowest survival limit. However, such a limit can also be reached by 
the absorption of J/ip by collisions with produced hadrons alone, as one can 
see from Fig. 4. It will be difficult to distinguish the two cases at such a high 
energy. It is therefore important to perform experiments at lower energies 
where such a distinction can be observed. 

5 Conclusions 

Using the relation between Et and b obtained from the Drell-Yan distribu- 
tion, we have analyzed the transverse-energy dependence of J /ip production. 
We found that the anomalous suppression of J/ip in Pb-Pb collisions can be 
explained by a model in which a new phase of strong absorption sets in when 
the local energy density exceeds 4.2 GeV/fm^. The critical energy density ecrit 
which corresponds to the onset of the new phase is close to the quark-gluon 
plasma energ3z_.density, ec ^ 4.17 GeV/fm'^, calculated from the lattice gauge 
theory result Ej of tc/T^ ~ 20 with Tc ^ 0.2 GeV. Therefore, it is interesting 
to speculate whether the new phase of strong absorption may be the quark- 
gluon plasma. However, other signatures are needed to corroborate whether 
the quark-gluon plasma has been produced. 

Our formulation of the absorption model in terms of produced hadron 
densities allows us to calculate J /ip survival probabilities at higher energies. 
We found interesting new results that the J /ip survival probability reaches the 
lowest survival limit when the collision energy ^/s reaches about 35 GeV. At 
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the energy of ^/s = 200 GeV, it will become difficult to distinguish absorption 
by the new phase or absorption by produced soft particles, as both effects will 
lead to the lowest survival limit. 
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